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OBSERVATION  OF  MAGNETIC  FIELDS  IN  LASER-PRODUCED 
PLASMA  USING  THE  ZEEMAN  EFFECT 


I.  INTRODUCTION 

It  has  been  known  for  some  time  that  a  large  magnetic  field1  is  generated  when  energetic  laser 
pulses  irradiate  solid  targets.  These  fields  can  result  from  thermoelectric  and  pressure  effects  that  are 
present  in  these  laser-generated  plasmas,  and  are  a  strong  function  of  temperature  and  pressure  gra¬ 
dients  within  the  plasma.  Knowledge  of  these  fields  is  important  in  studies  of  laser  fusion  and  inter¬ 
streaming  plasma  instabilities  due  to  their  effect  on  thermal  transport,2  absorption,  and  other  processes 
that  occur  in  these  plasmas.  For  example,  recently  Forslund  and  Brackbill1  have  predicted  that  these 
large  magnetic  fields  will  grossly  alter  the  transport  of  hot  electrons,  energy  and  ions  in  C02  laser-target 
interactions.  A  variety  of  techniques  to  measure  the  magnetic  fields  have  already  been  used,  such  as 
small  probe  coils,4  Faraday  rotation,*'6  and  the  imprint  on  magnetic  tape.7.  However,  this  paper 
describes  the  first  application  of  the  Zeeman  effect  to  measure  these  self-generated  magnetic  fields  in 
laser-plasmas,  and  the  first  observations  in  the  range  of  100-300  kG.  Those  fields  are  produced  by  laser 
irradiances  in  the  range  of  1012-1014  W/cm2  and  are  important  in  laser  ablation  studies. 

Although  magnetic  probe  coils  were  first  used  to  measure  magnetic  fields  in  laser-inadiated  solid 
targets,  they  can  only  be  deployed  at  distances  down  to  a  few  mm  from  the  target  surface  (due  to  the 
fact  that  the  plasma  severly  damages  the  probe)  and  consequently  could  only  measure  fields  of  a  few 
kilogauss.  Also,  the  magnetic  probes  perturb  the  plasma.  On  the  other  hand,  optical  measurements 
such  as  Faraday  rotation  and  Zeeman  effect  do  not  perturb  the  plasma  and  have  the  ability  to  measure 
the  magnetic  fields  at  much  higher  densities.  Stamper  et  al.5  and  Raven,  Willi  and  Rumsby6  have  used 
Faraday  rotation  to  measure  the  magnetic  fields  in  laser-produced  plasmas  at  laser  irradiances  of 
10l5-1016  W/cm2  and  have  reported  fields  of  1-2  megagauss  (MG).  In  those  measurements,  density 
gradients  limited  the  distance  that  the  probing  light  beam  could  penetrate  the  plasma  to  values  greater 
than  100  fxm  from  the  target  surface. 

A  comparison  of  the  advantages  and  disadvantages  of  Faraday  rotation  and  Zeeman  splitting  to 
measure  magnetic  fields  will  now  be  given.  Faraday  rotation  (1)  requires  an  external,  polarized  light 
source,  (2)  also  time  resolution  (typically  less  than  1  nsec)  if  the  external  light  source  is  of  short  dura¬ 
tion,  (3)  gives  information  that  is  the  integral  of  the  product  of  magnetic  field  and  electron  density 
along  the  optical  path  of  the  probe  beam,  and  (4)  gives  good  spatial  resolution  transverse  to  the  probe 
beam.  On  the  other  hand,  Zeeman  effect,  which  produces  the  splitting  and  broadening  of  spectral  lines 
emitted  by  the  plasma,  (1)  is  a  passive  diagnostic,  (2)  allows  time-resolved  measurements  (typically  a 
few  nsec)  to  be  made,  (3)  gives  the  average  magnetic  field  (but  not  the  electron  density)  along  the  line 
of  sight,  (4)  gives  moderate  spatial  resolution  of  the  magnetic  field,  and  (S)  with  the  proper  choice  of  a 
spectral  line  allows  observations  to  a  higher  density  region  of  the  plasma  than  would  the  Faraday  rota¬ 
tion  method.  Unfortunately,  the  Zeeman  splitting  technique  requires  considerable  effort  to  analyze  the 
spectral  line  shape,  since  other  broadening  mechanisms  such  as  thermal  Doppler  and  Stark  broadening 
as  well  as  instrumental  broadening  and  Doppler  shifts  due  to  mass  motion  are  also  present. 

Peacock  and  Norton*  had  used  the  Zeeman  effect  to  measure  megagauss  magnetic  fields  in  a 
plasma  focus  device.  For  these  measurements  they  chose  the  helium-like  CV  ls2s3Si-ls2p1/>2.i.o  raul- 
tiplet  spectral  lines  which  were  emitted  from  their  plasma.  We  have  chosen  the  same  muftiplet  with 
spectral  lines  for  J  -  2  at  2270.9  A,  J  -  1  at  2277.9  A,  and  J  -  0  at  2277.3  A  for  the  Zeeman  splitting 
measurements  described  in  this  report.  Reasons  for  the  choice  of  this  multiplet  are  (1)  it  is  relatively 
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insensitive  to  Stark  broadening,  since  the  line  is  a  transition  from  energy  states  with  low  principal  quan¬ 
tum  numbers,  (2)  observations  could  be  made  in  the  quartz  UV  where  standard  polarization  sensitive 
optics  can  be  used,  (3)  the  spectral  line  comes  from  a  relatively  high  temperature  (—100  eV)  region  of 
the  plasma  where  one  expects  the  magnetic  field  to  be  present,  and  (4)  calculation  of  Zeeman  splitting 
in  the  intermediate  magnetic  field  region  was  tractable.  LS  coupling  could  be  assumed  in  ouc  calcula¬ 
tions  of  the  transition  from  weak  field  (anomalous  Zeeman  effect)  to  strong  field  effect  (Paschen-Back9 
effect). 

Irons,  McWhirter  and  Peacock10  did  extensive  spectroscopic  investigations  of  the  C+4  ions  pro¬ 
duced  by  the  laser  irradiation  of  polyethylene  foils.  In  these  studies,  they  showed  that ..  j  C+4  ions  can 
be  seen  up  to  a  centimeter  from  the  target  surface,  and  further,  they  noted  that  although  the  ions  pri¬ 
marily  are  emitted  perpendicular  to  the  target  surface,  there  is  a  radial  velocity  component  (several 
times  smaller  than  the  axial  component)  parallel  to  the  target  surface.  Also,  Herbst  and  Grun11  in  a 
recent  study  using  a  novel  tracer-dot  technique  have  shown  that  the  fluid  flow  of  ablated  plasma  from 
planar  target  curves  away  from  the  target  normal  with  appreciable  radial  velocity  component. 

Further,  it  is  shown  in  this  report  that  when  absorption  of  the  spectral  line  occurs  in  conjunction 
with  the  Zeeman  effect,  it  is  possible  to  measure  smaller  magnetic  fields  at  a  given  wavelength  resolu¬ 
tion  than  could  be  measured  without  absorption.12  In  this  study,  magnetic  fields  were  measured  with 
field  strengths  of  200  kilogauss  (kG)  and  below. 


II.  SPECTRAL  LINE  BROADENING  EFFECTS 
A.  General  Considerations 

According  to  Zeeman  effect  theory,13  when  a  plasma  is  in  a  magnetic  field,  the  light  emitted  per¬ 
pendicular  to  the  magnetic  field  lines  is  linearly  polarized  with  one  set  of  components  polarized  parallel 
to  the  magnetic  field  (ir-components)  and  the  other  set  of  components  polarized  perpendicular  to  the 
magnetic  field  (<r -components).  In  this  experiment  observations  will  be  made  perpendicular  (as  near  as 
possible)  to  the  azimuthal  magnetic  field  lines;  by  using  polarizers  the  profiles  of  the  spectral  lines  con¬ 
taining  the  Jr -components  can  be  separated  from  those  containing  the  <r -components.  To  obtain  a 
measurement  of  the  magnetic  field  using  Zeeman  splitting  one  must  address  also  the  other  processes 
that  lead  to  the  broadening  or  wavelength  shifting  of  the  spectral  lines.  Generally,  it  is  necessary  to 
have  the  width  of  the  Zeeman  splitting  comparable  to  or  greater  than  the  net  broadening  from  these 
other  processes. 

Since  the  plasma  emitting  the  CV  multiplet  is  a  high  density,  hot  plasma  one  has  to  consider  the 
broadening  of  the  spectral  lines  by  Stark  and  thermal  Doppler  effects.  To  estimate  the  broadening  of 
lines  by  Stark  effect,  we  have  used  interferometric  measurements14  made  under  similar  conditions  to 
get  an  estimate  of  the  electron  density  and  then  used  the  following  relation  for  the  broadening  caused 
by  inelastic  and  superelastic  collisions;13 

BAj-dO*1)^  [y]  geNe.  (1) 

Here  #AS  is  the  full  width  at  half-intensity  in  A,  X  is  the  wavelength  of  the  line  in  A,  c  is  the  velocity 
of  light,  q,  is  the  excitation  rate  for  the  transition  ls2s3^  — *  ls2p3P  (Id-9  cm3s“’)  which  was  calculated 
by  Blaha,1*  and  N,  is  the  electron  density  in  cm-3.  The  excitation  rate  is  a  weak  function  of  tempera¬ 
ture  and  is  found  to  be  dominated  by  electron  collisions.  (Should  the  electron  temperature  be  lower  by 
a  factor  of  2,  the  Stark  width  would  increase  by  no  more  than  20  percent)  At  the  time  and  position  of 
observation  in  our  experiment  the  electron  density  is  estimated  to  be  between  lO^-lO20  cm'3.  Using 
these  values  in  Eq.  (1),  the  full  width  at  half-intensity  of  the  Stark  broadened  line  lies  between  0.18 
and  1.8  A  at  an  electron  temperature  of  100  eV. 
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The  thermal  Doppler  broadening  of  the  spectral  lines  can  be  estimated  if  the  ion  energies  are 
known.  The  electron  temperature  is  estimated  to  be  between  80-100  eV,  which  is  the  temperature 
range  for  maximum  emission  of  the  CV  lines,16  and  the  ion  temperature  should  approximate  the  elec¬ 
tron  temperature  when  the  electron  densities  are  this  high.  The  following  equation  can  then  be  used  to 
calculate  the  width  of  a  spectral  line  from  thermal  Doppler  broadening;17 

8Xj>  -  7.16  x  10-7X(r/A/)1/2,  (2) 

where  6X0  is  the  full  width  at  half-intensity  in  A,  T  is  the  ion  temperature  in  keV,  X  is  the  wavelength 
of  the  spectral  line  in  A,  and  M  is  the  atomic  weight.  For  an  estimated  ion  temperature  of  80-100  eV, 
the  thermal  width  of  the  CV  2270.9  A  line  would  fall  between  0.44  A  and  0.5  A.  Also,  if  the  instru¬ 
mental  width  is  known,  an  estimate  of  the  Stark-Doppler  width  of  the  line  can  be  made  from  the  n- 
profile,  which  is  only  weakly  sensitive  to  the  magnetic  field  strength.  This  estimate,  including  a  correc¬ 
tion  for  opacity  effects  (described  later  in  this  paper),  gave  a  Stark-Doppler  width  for  these  conditions 
of  0.6  A,  which  falls  between  the  values  derived  from  the  estimated  Nt  and  7}  given  above. 

In  addition  to  the  thermal  Doppler  width  of  the  spectral  line,  there  is  an  actual  shift  of  the 
wavelength  of  the  line  due  to  mass  motion.  This  is  due  to  the  radial  motion  of  the  plasma  plume  away 
from  the  laser  axis.  These  velocities  are  about  35%  of  the  axial  velocities.  Normally,  in  an  optically 
thin  plasma  this  would  give  a  wavelength  shift  both  to  the  red  and  to  the  blue  since  the  expansion 
occurs  both  toward  and  away  from  the  observer.  However,  as  will  be  explained  later  in  this  report,  only 
the  blue-shifted  emission  is  seen  because  the  red-shifted  emission  has  to  pass  through  the  dense  central 
region  of  the  plasma  plume  and  is  almost  completely  absorbed.  The  wavelength  shift  was  measured  to 
be  1 .4  A  to  the  blue.  Using  the  Doppler  formula, 

-y-  -v/c.  (3) 

we  find  that  the  mean  radial  velocity  is  1.4-1.8  x  107  cm/sec,  which  is  consistent  with  other  measure¬ 
ments.11  This  causes  little  problem  in  the  interpretation  of  the  data,  because  all  of  the  lines  measured 
have  a  constant  wavelength  shift,  which  can  simply  be  subtracted  from  the  data. 

The  instrumental  width  of  the  spectrograph  was  determined  by  scanning  a  very  narrow  spectral 
line  from  a  Geissler-tube  light  source.  It  was  found  that  the  spectral  line  shape  was  a  Gaussian  with  the 
full  width  at  half-intensity  corresponding  to  the  widths  of  the  entrance  and  exit  slits  of  the  monochro¬ 
mator.  Two  slit  widths,  200  and  100  pm,  were  used  in  this  study,  giving  instrumental  widths  of  1.5  A 
and  0.75  A,  respectively.  Having  some  confidence  in  the  shape  of  the  broadened  spectral  lines,  it  is 
now  possible  to  use  this  data  to  calculate  the  shapes  of  the  multipiet  for  different  values  of  the  mag¬ 
netic  field. 

First,  it  is  necessary  to  calculate  positions  and  intensities  of  the  ir  and  <r  components  of  the  CV 
multipiet  for  the  range  of  magnetic  fields  that  we  are  considering,  i.e.,  0-1  MG.  Since  we  are  in  the 
intermediate  magnetic  field  strength  regime,  the  weak-field  approximation  cannot  be  used  and  a  more 
detailed  calculation  must  be  done.  This  calculation  is  described  in  the  Appendix.  A  typical  Zeeman 
pattern  for  B  —  0.5  MG  is  shown  in  Fig.  1.  The  scale  along  the  abscissa  is  in  wave  numbers  and  0  is 
located  at  the  center  of  the  /  -  2  (7S(  -  3P7)  line  (2270.9  A).  The  v  components  are  plotted  above 
the  abscissa  and  the  <r  components  are  plotted  below  the  abscissa.  Both  the  position  and  the  relative 
intensity  of  each  component  can  be  determined  from  this  plot.  Immediately,  it  is  noted  that  there  is 
appreciable  asymmetry  compared  with  the  weak-field  anomalous  Zeeman  splitting  patterns,  e.g.,  note 
the  difference  in  the  w  intensities  of  the  M  —  1  and  Af  -  — 1  components  on  the  two  sides  of  the  /  —  2 
line. 
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From  these  calculated  Zeeman  intensity  patterns  of  the  components,  a  calculation  can  be  done 
which  will  give  the  expected  shape  of  the  multipiet  as  viewed  in  the  two  polarizations,  *  and  <r.  (See 
Fig.  2.)  To  do  this  a  small  computer  code  was  written  in  which  each  individual  component  of  the  Zee- 
man  pattern  is  broadened  by  the  net  amount  of  Stark,  thermal  Doppler  and  instrumental  broadening 
(but  without  opacity  effects),  and  then  the  individual  contributions  of  all  components  are  added 
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j  -ZEEMAN  COMPONENT  WITH  B*0.5MG 
Fig.  1  —  Zeeman  pattern  for  the  triplet  CV  JS,  -  ’/‘u  #  at  B  -  0.5  MO  ($0  D 

together  to  give  the  emitted  optically-thin  line  spectral  shape.  These  computed  lines  are  then  similar  in 
shape  to  the  spectral  lines  which  are  actually  measured  and  by  choosing  the  calculated  profile  which  best 
fits  the  measured  profile  one  can  determine  the  approximate  value  of  the  magnetic  field  strength. 

B.  Opacity  Corrections 

The  observed  shapes  and  ratios  of  the  <r  to  w  line  intensities  suggested  that  opacity  effects  were 
important.  It  is  well  known  that  the  absorption  coefficient  of  a  spectral  line  is  strongly  dependent  on 
the  line  shape,19  and  at  zero  magnetic  field  is  a  maximum  at  the  center  of  the  spectral  line.  Thus,  the 
absorption  coefficient  at  the  line  center  will  be  different  for  the  tr  and  a  profiles  in  the  presence  of  a 
magnetic  field.  This  is  because  the  ir  components  peak  at  the  same  wavelength  as  the  spectral  line  at 
zero  magnetic  field,  whereas  the  <r  components  have  a  dip  at  the  center  of  the  line.  When  the  opacity 
of  the  plasma  is  appreciable,  the  n  components  will  therefore  suffer  more  absorption  than  the  <r  com¬ 
ponents.  This  will  change  the  total  intensity  of  the  line  so  that  the  magnetic  field  effects  become 
observable  at  relatively  low  fields,  i.e.,  even  when  differences  in  the  optically  thin  line  shapes  are  still 
obscured  by  instrumental  broadening.  Indeed,  the  relative  amplitudes  of  the  jt  and  <r  profiles  can  be 
inverted. 

Opacity  effects  in  the  experiment  within  the  annular  emitting  region  were  modelled  by  setting  the 
opacity  of  the  emitting  region  ( kKl )  equal  to  C|/,(X),  where  I,(\)  is  the  total  thin-line  intensity  or  line 
shape.  Opacity  effects  in  the  cooler,  outer  region  were  similarly  modelled  by  setting  the  opacity  of  the 
outer  region  equal  to  C2/,(X).  The  coefficients  Ct  and  C2  are  adjustable  parameters.  Since,  in  the 
code,  the  total  intensity  7,(X)  of  all  components  is  separately  normalized  to  unity  for  the  w  and  a  com¬ 
ponents,  opacity  effects  will  only  become  dominant  when  the  coefficients  are  much  greater  than  one. 
The  relation  giving  the  emitted  intensity  of  the  line  with  opacity  effects  included  is, 

7*(x)  -  jr  |l  -  exp  l-C,7,(X)j]  exp  -  fC2/,(x)J.  (4) 

c>  I  Emission  Region  [  Outside  Region 

It  was  found  in  matching  the  data,  that  the  coefficient  C i  for  the  emitting  region  was  not  large  com¬ 
pared  to  one.  Also,  the  opacity  correction  factor  1  -  exp  t-C|/,(x)l  for  the  emitting  region  was  nor¬ 
malized  (by  dividing  by  C()  so  that  the  optically  thin  intensity  7,(x)  is  recovered  in  the  limit  of  small 
C|.  However,  a  large  value  of  C2  (—5)  was  required  to  match  the  data;  the  correction  factor 
exp  (~C2/,(X)1  is  the  fraction  of  light  transmitted  through  the  cooler  absorbing  (outer)  region.  Note 
that  when  both  regions  are  optically  thin  (both  Ct,  C2  small),  the  optically  thin  emission  /,(X)  is 
recovered.  The  opacity  effect  is  included  in  the  analysis  before  the  instrument  function  is  convolved 
with  7,(x)  to  get  /,'(X),  which  can  then  be  directly  compared  with  the  observed  line  profiles. 
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Fif.  2  —  Calculated  broadened  tinea  including  Zeeman 
ipUtting  of  the  CV  multiple!  ai  a  function  of  magnetic  field 
strength  in  the  range  0-1000  kO.  The  w  profile  Unet  are 
plotted  above  the  abectaaa  and  the  <r  profile  below  the 
abscissa.  Thermal  and  Stark  broadening  (1.6  A)  and  in¬ 
strument  broadening  (0.7$  A)  are  included  in  the  calcula¬ 
tion. 
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III.  EXPERIMENT 


In  this  experiment,  a  pulsed  Nd:phosphate  glass  laser  beam  (30-40  J  in  5  nsec)  at  1.0S  ftm  irradi¬ 
ates  a  50- ftm  thick  carbon  target  in  the  near  held  of  an  f/6  lens.  (See  Fig.  3a.)  In  the  first  set  of  the 
measurements  a  7-cm  diameter  mask  was  placed  ahead  of  the  focusing  lens  so  that  a  ring-shaped  laser 
intensity  pattern  with  an  outer  diameter  of  approximately  0.3S  mm  and  a  deep  minimum  would  fall 
onto  the  target  (allowing  a  larger  density  gradient  and  magnetic  field  to  develop).  In  the  second  set  of 
measurements  the  mask  was  removed  from  the  beam  so  that  a  focal-intensity  distribution  which  is 
almost  flat-topped  (but  with  a  shallow  minimum  at  the  focal  spot  center)  would  be  incident  on  the  tar¬ 
get.  The  plasma  light  was  collected  by  an  f/ 1.6  parabolic  mirror.  Ml.  The  mirror  was  placed  one  focal 
length  from  the  target  so  that  a  collimated  light  beam  would  be  transmitted  via  plane  mirrors,  M2,  M3, 
M4  to  the  f/S  concave  mirror  MS  which  focused  the  beam  through  a  Wollaston  prism  onto  the  slit  of  a 
1-m  spectrograph.  The  collected  light  always  had  small  angles  of  incidence  with  respect  to  the  mirror 
normals  in  order  to  minimize  depolarization  effects.  The  magnification  of  this  optical  train  was  S.3  at 
the  entrance  slit  of  the  spectrograph.  The  optics  were  aligned  so  that  the  image  of  the  spectrograph  slit 
at  the  center  of  the  plasma  plume  was  parallel  to  the  target  surface  and  400  ^m  from  the  surface.  (See 
Fig.  3b.)  (This  distance  was  chosen  because  a  reasonable  spectral  line-to-continuum  ratio  could  be 
obtained  there.  At  distances  closer  to  the  target  surface,  the  continuum  intensity  overpowered  the  line 
intensities.)  The  cross-sectional  area  of  the  slit  image  at  the  center  of  the  plume  is  —20  ^m  (or  40  ;j.m 
for  the  200  ftm  slit  width)  by  —200  fim,  and  emitted  light  from  f/1.6  cones  on  both  sides  of  the  focal 
region  is  accepted  by  the  detector. 

The  Wollaston  prism,  shown  in  the  side  view  of  Fig.  3a,  vertically  separated  (along  the  slit)  the 
light  beam  from  the  plasma  plume  into  the  two  orthogonal  polarizations,  n  and  «r,  so  that  these  polari¬ 
zations  would  then  continue  to  be  separated  in  passing  through  the  spectrograph.  Mirrors  at  the  exit 
slit  of  the  spectrograph  then  directed  the  two  polarized  signals  into  separate  photomultiplier  tubes.  To 
optimize  the  intensity  collected  by  the  detectors  at  this  wavelength,  a  spectrograph  grating  blazed  at 
2000  A  was  used.  The  signal  from  the  photomultipliers  (RCA  1P28)  was  sent  through  50-0  cables  to 
fast  recording  oscilloscopes  (Tektronix  7104),  giving  a  system  time-response  of  3  nsec.  A  relative 
intensity  calibration  of  the  spectrograph  and  detectors  was  made  using  a  deuterium  arc  lamp.  Also,  the 
continuum  intensity  at  2256  A,  which  should  not  be  polarized,  was  used  to  calibrate  on  a  relative  scale 
the  entire  optical  train.  A  calibration  of  the  spectrograph  wavelength  scale  was  made  using  a  Hg  Pen- 
Ray  light  source.  Using  this  optical  setup  the  intensities  of  the  it  and  a  emissions  at  a  particular 
wavelength  are  recorded  separately  and  simultaneously.  To  obtain  the  spectral  line  shapes  it  is  neces¬ 
sary  to  assume  shot-to-shot  reproducibility  and  to  take  a  sufficient  number  of  shots  to  cover  the  line 
shape,  one  smalt  increment  of  wavelength  at  a  time. 
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(b) 

Fig.  3  —  (a)  Experimental  arrangement  (b)  Area  of  observation 


IV.  RESULTS  AND  ANALYSIS 
A.  General 

Typical  photomultiplier  oscillograms  are  shown  in  Fig.  4.  The  two  traces  shown  were  taken  on  the 
same  shot  near  the  peak  of  one  of  the  CV  lines;  (a)  signal  with  7-polarization  and  (b)  signal  with  cr- 
polarization.  Two  prominent  peaks  are  seen  on  each  trace.  The  first  peak  is  due  to  the  early  rise  of  the 
continuum  (and  possibly  the  CV  line)  intensity.  Since  the  peak  temperature  that  is  produced  in  the 
plasma  plume  is  greater  than  100  eV,  e.g.,  about  200  eV,  the  CV  line  initially  bums  through  and  does 
not  radiate  significantly  until  the  temperature  falls  later  in  time;  then  recombination  occurs  and.  the  CV 
line  reappears.  Thus,  the  second  peak  is  primarily  due  to  the  CV  line  emission  as  the  plasma  tempera¬ 
ture  decreases  to  between  80  and  100  eV.  (The  contribution  of  the  underlying  continuum  measured  on 
a  separate  shot,  is  subtracted  off  of  the  signal  shown  here.)  A  fiducial  signal  is  fed  into  the  oscilloscope 
so  that  the  time  of  the  CV  signal  with  respect  to  the  time  of  the  peak  of  the  laser  pulse  can  be  deter¬ 
mined.  Near  the  spectral  line  intensity  peak,  the  line-to-continuum  ratio  is  3  or  more. 

It  is  assumed  in  these  studies  that  the  CV  (C**)  ions  are  distributed  in  a  conical  annulus  about 
the  target  normal.  This  is  expected  for  the  case  of  a  ring-shaped  focal  spot,  but  is  also  produced  from 
the  intensity  distribution  with  no  mask  in  the  beam  as  well;  due  to  our  intensity  distribution  having  a 
shallow  central  minimum.  Figure  5  shows  a  cross  sectional  view  of  the  plasma  plume.  In  this  r  the 
interior  of  the  plume  is  cooler  and  more  dense  than  the  surrounding  ring  (due  to  the  lower  lasr  ten¬ 
sity  at  the  center  and  the  flow  of  mass  from  the  ring  into  the  central  region.)1114  Because  the  'ilar 


Fig.  4  —  Typical  photomultiplier  oscillograms.  The  time  r  -  0 
corresponds  to  the  time  of  the  peak  of  the  laser  pulse. 
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Fig.  S  —  Cross-sectional  view  of  plasma  plume 


ring  is  expanding  away  from  the  center,  an  observer  would  see  photons  produced  on  the  near  side  of 
the  ring,  blue-shifted,  and  those  on  the  far  side  of  the  ring,  red-shifted.  However,  because  the  center 
region  of  the  plume  is  relatively  cool,  is  at  a  higher  density  and  is  large  (~mm),  the  red-shifted  pho¬ 
tons  are  absorbed  and  do  not  reach  the  observer.  As  discussed  earlier,  the  blue-shifted  CV  lines  are 
also  partially  absorbed,  particularly  at  later  times,  by  the  region  outside  of  the  annular  region  emitting 
the  CV  light. 

Although  the  magnetic  field  is  in  a  clockwise  direction  in  the  drawing  in  Fig.  5,  the  Zeeman  effect 
cannot  distinguish  a  clockwise  from  a  counterclockwise  field.  Observations  of  the  line  emission  are 
made  perpendicular  to  the  laser  axis  and  intersecting  the  laser  axis,  so  that  the  observations  are  approxi¬ 
mately  perpendicular  to  the  magnetic  field. 

As  mentioned  in  the  introduction,  to  determine  a  magnetic  field  from  the  line  profile  measure¬ 
ments,  one  must  assume  the  composite  shape  of  the  Stark  and  Doppler  broadening  of  the  spectral  line, 
use  the  measured  instrumental  width,  and  calculate  the  Zeeman  pattern  for  the  expected  range  of  the 
magnetic  field.  If  absorption  is  also  present  (as  is  usually  the  case),  this  too,  will  have  to  be  taken  into 
account.  Then  one  can  determine  the  magnetic  field  which  will  give  the  best  fit  to  the  lin  profile. 
Since  this  analysis  involves  a  rather  tedious  calculation,  an  interactive  computer  code  was  written  for 
this  purpose.  Our  calculation  utilized  the  values  of  shifts  and  intensities  for  intermediate  magnetic 
fields  of  each  ir  and  <r  component  as  described  in  the  Appendix.  The  steps  in  the  calculation  are  as  fol¬ 
lows: 

1.  First,  a  Stark-Doppler  profile  of  the  individual  Zeeman  components  was  assumed.  (This  decision 
was  guided  by  estimates  of  Ne  from  interferometric  measurements  and  T,  from  the  CV  emission  data.) 

2.  This  Stark-Doppler  width  was  then  assigned  to  both  the  it  and  <r  components  of  the  Zeeman  profile 
for  a  particular  assumed  value  of  the  magnetic  field.  Each  component  is  shifted  from  its  zero  field  loca¬ 
tion  and  its  relative  amplitude  set  according  to  its  intermediate  field  calculation.  The  profiles  of  the  3  ir 
components  of  the  J  -  2  line  (or  the  6  <r  components)  are  then  added  together  to  give  a  net  ir  (or  or 
profile);  at  this  point  no  absorption  effects  have  been  included.  We  now  have  the  optically  thin  profile 
/,(A). 

3.  Next,  the  opacity  effects  are  included.  This  is  done  for  the  emitting  and  outside  regions  of  the 
plasma  by  adjusting  coefficients  Ct  and  C2  and  using  Eq.  (4). 
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4.  Finally,  then,  to  calculate  the  "observed”  line  profiles,  one  must  convolve  the  /,(x)  profile  with  the 
known  instrument  profile.  The  instrument  function  used  was  a  Gaussian  with  a  FWHM  of  l.S  A,  or 
0.7 5  A,  as  appropriate. 

5.  The  computed  profile  /«  (X)  now  can  be  compared  to  the  measured  profile;  if  it  does  not  fit,  a  new 
value  of  the  magnetic  field  can  be  tried,  a  new  Stark-Doppier  width  can  be  used,  or  the  values  of  the 
opacity  C  i  and  C2  may  be  varied  to  obtain  a  satisfactory  fit. 

B.  Experiments  with  a  Ring-Shaped  Laser  Focal  Spot 

In  Fig.  6  a  comparison  is  shown  of  the  experimental  profile  for  the  /  -  2  line  of  the  CV  multiplet 
with  a  theoretical  profile  for  200  kilogauss.  Since  the  data  for  J  -  0,  J  -  1  lines  is  incomplete  we  did 
not  include  them  in  the  paper.  The  experimental  profile  was  observed  S  nanoseconds  after  the  peak  of 
the  laser  pulse,  and  near  the  time  of  the  peak  intensity  of  the  J  -  2  line  (centered  at  2270.9  A).  The 
measured  data  points  of  the  w  components  are  noted  with  crosses  and  the  data  points  for  the  a  com¬ 
ponents  are  shown  with  circles.  We  assume  reproducibility  of  the  plasma  plume  and  take  only  one  or 
two  shots  for  each  wavelength  (which  gives  simultaneous  data  for  both  it  and  <r  components).  The 
scatter  in  the  data  is  indicative  of  the  error  in  this  procedure.  Generally,  the  data  points  fall  within  a 
scatter  bracket  of  ±  20%;  however,  at  +  3  A  the  scatter  is  higher.  (This  is  possibly  due  to  the  fact  that 
the  /  —  0  line  is  beginning  to  overlap  the  J  —  2  line  at  that  wavelength.)  For  the  fit  shown  here,  rather 
small  opacity  parameters  C\  -  0.3  and  C2  -  2  are  required.  The  agreement  in  profile  fits  was  better  at 
200  kG  than  at  100  kG  or  300  kG,  although  it  was  noted  that  a  slightly  smaller  n  width  would  agree 
better  with  the  data,  indicating  that  a  field  slightly  larger  than  200  kG  would  agree  somewhat  better 
with  the  data. 


AT  t  =.  +5  nsec  THEORETICAL  CURVES 


X  (ANGSTROMS) 

Fig.  6  —  Comparison  of  experimental  data  and  calculated  profiles  for  J  -  2  line  at  t  —  5 
nsec  for  ring-shaped  laser  focal  spot.  Measured  data  points:  X  are  for  w  polarization  and 
0  are  for  a-  polarization.  Curves  are  calculated  using  an  instrument  width  “  1.5  A,  the 
Stark-Doppier  width  —  1.44  A,  and  B  -  200  kG. 
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Data  was  also  analyzed  at  a  later  time  (8  nsec  after  the  peak  of  the  laser  pulse)  on  these  same 
shots.  The  profile  obtained  at  this  time  is  show>.  .r  Fig.  7a.  The  striking  feature  in  this  data  is  that  the 
<t  line  is  nearly  twice  as  intense  as  the  ir  line  at  the  center  of  the  profiles.  It  was  found  that  by  making 
the  surrounding  region  more  absorbing  (C2  -  5)  and  continuing  to  assume  a  magnetic  field  of  200  kG 
the  theoretical  curve  gives  a  fairly  good  fit  to  the  experimental  data.  A  slightly  smaller  or  width  would 
give  an  even  better  fit  to  the  data,  indicating  that  the  field  is  just  under  200  kG.  The  important-  point  to 
be  made  regarding  this  data  is  that  the  only  way  the  o-  profile  could  have  a  much  higher  intensity  than 
the  ir  profile  is  for  absorption  to  be  taking  place,  in  which  the  cr  profile  in  the  presence  of  a  magnetic 
field  has  a  dip  in  the  line  center  (where  peak  absorption  would  normally  occur)  and  therefore,  has  less 
absorption  than  the  ir-profile  which  has  a  sharp  peak  at  the  line  center  even  with  a  magnetic  field. 

1.  Error  Analysis 

We  have  several  free  parameters  which  can  be  varied  to  get  a  fit  to  the  experimental  data.  Since 
we  do  not  have  precise  auxiliary  data  to  pin  down  these  variables,  it  is  necessary  to  do  trial  calculations 
varying  the  parameters  individually  to  see  if  we  can  get  a  satisfactory  fit  with  different  sets  of  values  of 
the  parameters,  i.e.,  to  see  if  we  have  a  unique  fit  with  the  values  chosen.  For  this  set  of  calculations, 
the  value  of  B  was  varied  between  0  and  300  kG,  Cx  was  varied  from  0.3  to  7,  C2  was  varied  from  0.1 
to  10,  and  the  Stark-Doppler  width  (WS.D)  was  varied  from  0.36  A  to  1.5  A  (where  the  upper  limit 
could  be  determined  by  the  width  of  the  ir  component,  whose  width  is  much  less  dependent  on  the 
magnetic  field  than  the  cr  component). 

From  over  40  calculations,  we  found  that  the  best  fit  occurred  with  Ws-D  —  0.6  A,  C]  —  0.3, 
Ci  —  5,  and  the  magnetic  field  —  200  kG.  An  example  of  the  sensitivity  of  the  line  shape  to  the  mag¬ 
netic  field,  holding  the  other  variables  constant  to  the  values  given  above  is  shown  in  Fig.  7b  and  7c. 
(All  of  the  curves  are  normalized  with  respect  to  the  a  values  at  0  A.)  It  is  noted  in  Fig.  7b  where 
B  —  0  G  that  the  it  and  o-  curves  are  identical,  and  thus  could  not  fit  the  data.  The  data  points  are  also 
plotted  on  the  same  set  of  curves  in  Fig.  7c.  It  is  seen  that  the  best  overall  fit  of  the  calculated  curves 
with  the  data  points,  considering  both  the  amplitude  and  the  shape  is  given  by  the  200  kG  <r  and  ir 
curves.  (A  somewhat  better  fit  might  be  obtained  if  increments  of  magnetic  field  smaller  than  100  kG 
were  used.)  From  this  error  analysis,  it  is  concluded  that  the  accuracy  for  the  magnetic  field  measure¬ 
ment  is  ±75  kG  when  absorption  is  present. 

C.  Experiments  with  a  Flat-topped  Focal  Spot  (Mask  Removed  from  Laser  Beam) 

When  the  mask  in  front  of  the  lens  is  removed,  the  focal  intensity  distribution  is  relatively  flat- 
topped  with  a  shallow  minimum  at  the  center.  The  diameter  of  the  focal  spot  was  0.48  mm,  which 
contained  90%  of  absorbed  laser  energy.  The  resulting  CV  data  is  shown  in  Fig.  8.  The  instrumental 
width  remained  1.5  A,  the  same  as  for  the  previous  data.  The  theoretical  curve  for  B  -  100  kG, 
Ci  -  .01,  C2  -  3.5  and  the  Stark-Doppler  width  of  1.5  A  fits  the  experimental  data  for  the  J  -  2  line 
fairly  well  in  both  shape  and  amplitude  for  both  the  ir  and  the  o-  lines.  Due  to  the  lack  of  structure  in 
the  <r  line,  it  appeared  desirable  to  have  better  wavelength  resolution.  For  this  reason  another  data  run 
was  taken  in  which  the  instrumental  width  was  decreased  from  1.5  A  to  0.75  A,  the  other  experimental 
conditions  remaining  approximately  the  same. 

The  results  with  the  improved  wavelength  resolution  are  shown  in  Fig.  9;  more  structure  can  be 
seen  in  the  profiles.  For  the  J  -  2  line  both  the  ir  and  the  <r  lines  are  in  good  agreement,  in  shape  and 
amplitude,  with  the  calculated  curve.  (The  Stark-Doppler  broadening  parameter  was  changed  from  1.5 
A  to  1.0  A  and  the  opacity  parameters  were  changed  from  Cj  -  .01  to  0.3  and  C2  -  3.5  to  3.0  from 
those  used  in  Fig.  8  since  the  wavelength  resolution  is  better  here  and  more  accurate  values  of  the 
above  parameters  can  be  chosen.)  Again  these  curves  have  been  plotted  for  different  values  for  B  to 
show  the  sensitivity  of  the  technique.  The  best  fit  was  again  B  -  100+/&  kG. 
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WAVELENGTH  (A) 

Fi|,  7  —  Comparison  of  measured  data  and  calculated  profile  for  7  —  2  line  at  r  —  I  naec 
for  the  same  data  runs  as  used  for  Fig.  6.  Measure  data  points:  X  are  for  *  polarisation 
and  0  are  for  tr  polarisation.  Curves  are  calculated  uaiug  an  instrumental  width  -  1.5  A, 
the  Stark -Doppler  width  -  0.6  A  and  the  B  varied  between  0  and  300  kO.  (a)  B  -  200 
kO,  (b)  B  -  0  O,  (c)  B  varied  between  100  kO  and  300  kO  to  show  error  spread  in  data. 
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V.  DISCUSSION  AND  CONCLUSIONS 

The  magnetic  field  produced  when  a  pulsed  laser  beam  is  focused  onto  a  target  is  measured  for 
the  first  time  using  the  Zeeman  effect.  These  measurements,  which  were  made  at  electron  densities  of 
1019  -  1020  cm-3,  allowed  somewhat  lower  magnetic  fields  (—100  kO)  to  be  measured  than  the  Faraday 
Rotation  technique,  and  further  did  not  require  the  measurement  of  the  electron  density. 

The  plasma  opacity  present  enhanced  the  accuracy  of  the  .0-field  measurement.  At  laser  irradi- 
ance  of  — 1012  W/cm2,  the  magnetic  field  increased  from  100  kO  to  200  kO  when  the  focal  spot  was 
changed  from  a  flat-topped  distribution  to  a  ring  pattern,  which  caused  higher  density  gradients.  The 
magnetic  field  measured  using  the  Zeeman  effect  with  the  ring  shaped  focal  spot  is  consistent  with  a 
Faraday  Rotation  measurement  made  under  similar  conditions. 

In  summary,  magnetic  fields  of  1-2  MG  were  seen  previously  when  the  laser  irradiance  was  high 
(— 1016  W/cm2)  and  pulse  duration  short  (100  psec);  here  the  magnetic  field  was  found  to  be  reduced 
to  200  kG  and  below  for  lower  irradiances  (S  x  1012  W/cm2),  longer  pulse  durations  (3-4  nsec)  and 
later  observation  times  (5-8  nsec). 
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Appendix 

THEORY  AND  CALCULATION  OF  ZEEMAN  EFFECT 
FOR  INTERMEDIATE  FIELDS 


The  field-free  ( ls2p)3Pj  levels  of  CV  for  J  -  0,  1  and  2  are  rather  closely  spaced,  with  intervals20 
Ej( 0)  -  £i(0)  —  136  cm-1  and  £](0)  -  £0(0)  -  -13  cm-1.  Because  of  the  small  (and  inverted)  inter¬ 
val  between  the  y  -  1  and  0  levels,  use  of  lowest  order  perturbation  theory  for  the  Zeeman  effect  as  in 
Ref.  8  would  not  be  justified  for  the  field-strengths  of  interest.  This  is  easily  seen  by  comparing  the 
magnitude  of  magnetic  interactions,21 

l^mj  =  H„B  “  46.686  B  cm-1,  (Al) 

with  B  in  megagauss,  with  the  above  field-free  intervals.  Particularly  important21  here  are  the  changes 
in  relative  intensities,  relative  to  the  low  field  (anomalous  Zeeman  effect)  results,  of  the  various  Zee- 
man  components  associated  with  the  2 iPJ  -  2 1Sl  transitions  observed  in  our  experiment. 


To  obtain  a  quantitative  theory  of  the  Zeeman  effect  for  intermediate  fields,  one  writes  the  2 iP 
eigen  functions  as  superpositions  of  LS eigen  functions  of  various  J  >  \M\, 

| LSM>  -  lCJit\LSJM> 

J  >  \M\  (A2) 

Only  the  magnetic  quantum  numbers  M  remain  'good,*  and  the  expansion  coefficients  Cju  follow  from 
the  Schrodinger  equations 

[fl}(0)  +  Haut]\LSM>  -  Em(B)\LSM>  (A3a) 

where  Hj(i 0)  is  the  field-free  Hamiltonian  with  eigenvalues  £}(0)  and  eigenfunctions  \LSJM> .  The 
nonvanishing  matrix  elements  of  the  magnetic  interaction  are  in  the  \LSJM>  representation22 

<  1 1  1 1  JM>  -  y  Mijl.B,  (A4a) 

<\  \  (/  —  1)M>  —  — y  ft0B  | (  — 1  (A4b) 


for  diagonal  and  off-diagonal  elements,  respectively.  According  to  the  various  properties  of  Hj( 0), 
/fmM,  and  I LSJM > ,  the  system  of  equations  corresponding  to  Eq.  (A3a)  becomes  for  each  M  a  set  of 
coupled  linear  equations 

£  (l£r(0)  -  Ery(B)]6rj  +  <y"3/|^m<#|yAf>}C/^(ff)  -  0.  (A3b) 


which  can  only  be  solved  if  the  level  energy  Ej<m(B)  is  such  that  the  determinant  of  the  matrix  in  curly 
brackets  vanishes.  For  each  such  root,  J',  one  then  finds  the  coefficient  vector  Cfy(B)  which  will  be 
assumed  to  be  normalized.  The  actual  matrices  are 


(£2(0)  -  Wff)  1  — 7?  HoB 

1 


0 


V5 

(£,(0)  -  £>«(£)) 

W2 

3 


HoB 


(£2(0)  ±  y  m.R  -  Eru(B)] 
-y  HoB 

l£2(0)  ±  3 fi»B  ~  Eru) 


[£„(0)  -  £>*(£)] 


-JH.B 


l£,(0)  ±  fn,B-  Erti(B)i 


(AS) 
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Fi*.  10  —  Zeeman  ifaifti  of  (J'.hf)  magnetic  aub-leveta  of  2>PJ  in  unite  of  f  «  1  and 
magnetic  perturbation  energy  n,M  -  46.616  cm-1  AMO).  To  obtain  the  actual  level 
positions.  e.g..  relative  to  the  zero- Ac  id  ]'  —  2  level,  123  cm-1  and  136  cm-1  must  be 
subtracted  from  the  /'  —  0  and  /  —  I  Zeeman  shifts.  This  procedure  causes  the  (0,0) 
level  to  stay  below  the  (2,0)  level  for  all  field  values,  i.e.,  there  is  no  crossing  of  levels. 


tor  Af  -  0,  hi  -  ±  1  and  hi  -  ±2,  respectively,  written  in  the  order  of  decreasing  indices  /and  /".  It 
is  convenient  to  refer  all  energies  in  the  numerical  calculations  to  the  unshifted  position  of  the  /  —  2 
level,  i.e.,  to  use  £2(0)  -  0,  £j(0)  —  -136  cm-1  and  £0(0)  —  -123  cm-3.  The  value  of  is  as  in 
Eg.  (Al). 

In  Fig.  10,  the  differences  of  the  eigenvalues  £>*(£)  and  EJtt(0)  are  plotted,  in  units  of  n,B, 
and  it  is  seen  that  significant  deviations  from  the  weak  field  limit  occur  for  B  >  0.1  MG.  (Note  that 
there  is  no  real  level  crossing  of  the  (0,0)  and  (2,0)  levels  near  S  MG.)  The  relative  intensities  of  the 
Zeeman  components  for  observation  perpendicular  to  the  field  can  be  calculated  from21 

A#  -  hi:  Si  -  Y.  i  hi  0  i|e<w  +  »i“i!  1 1  <AM 

V  -  M  ±  i  s;  -  i  1  T!  u  ±  ,Jl3(2Z  +  .)!«{;  ;  .  (A6b> 

where  the  expressions  in  large  parentheses  and  large  curly  brackets  are  Wigner’s  3  —  j  and  6-  j  sym¬ 
bols,  respectively.  The  6  -  J  symbols  have  the  values  1/3,  -1/3  and  1/3  for  /  -  0,  1  and  2.  They  may 
therefore  be  combined  with  the  phase  factors  into  a  common  factor  1/3.  Using  tabulated21  3  -  j  sym¬ 
bols  one  obtains  for  w- polarization 

J5  i 

hi -hi-  0;  Si  -  y-  Cfc  -  -j-  Cio  (A7a) 

hi  -  hi' -  ±\:  Si  -  j\Ci±{  *  Cf±,  |2 
hi-  hi  -  ±2:  Si  -0 
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(A7b) 

(A7c) 


Fig.  11  —  Relative  intensities  of  w  (up)  and  <r  (down)  Zeeman  components 
(J'.AfM)  of  the  23Pjl  —  23S,  transitions  u  functions  of  magnetic  field 
strength.  Observation  is  assumed  to  be  perpendicular  to  the  magnetic  field. 


and  for  tr  polarization 


0,  M'-  ±1 :  #  -  -jj  35  ^  cfo  +  c&l 

(A8a) 

M  -  ±  1,  M'  -  0:  Si :  -  ~\Ci±i  ±  Ci±x  I2 

(A8b) 

A/  -  ±  2,  ±\:SC  -  j-ICta!1-  j. 

(A8c) 

where  M'  is  the  magnetic  quantum  number  of  the  23S]  (lower)  state.  The  calculated  intensities  are 
shown  in  Fig.  11.  They  clearly  deviate  from  the  weak  field  limit  already  at  fields  below  0.1  MG. 

To  compute  line  rather  than  level  shifts,  the  shifts  of  the  lower  level,  namely 

Em ■  -  2 M'fi.B  -  93.372  M'B  cm'1  (A9) 

must  be  subtracted  and  the  zero-field  23P2  -  23Si  interval  be  added.  Or,  if  wavelengths  are  desired, 
one  calculates  the  wavelength  shifts  according  to 

AXjC* - “  MO)  -  Em  (B)\  (A10) 

with  X2  -  2270.9,  X,  -  2277.9  and  X0  -  2277.3  A. 
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HANSCOM  AFB,  MA  01731 

OICY  ATTN  EEP  A.  LORENTZEN 

DEPARTMENT  OF  ENERGY 
LIBRARY  ROOM  G-042 
WASHINGTON,  D.C.  20545 

OICY  ATTN  DOC  CON  FOR  A.  LABOWITZ 


DEPARTMENT  OF  ENERGY 
ALBUQUERQUE  OPERATIONS  OFFICE 
P.O.  BOX  5400 
ALBUQUERQUE,  NM  87115 

OICY  ATTN  DOC  CON  FOR  D.  SHERWOOD 

EG&G,  INC. 

LOS  ALAMOS  DIVISION 

P.O.  BOX  809 

LOS  ALAMOS,  NM  85544 

OICY  ATTN  DOC  CON  FOR  J.  BREEDLOVE 

UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.O.  BOX  808 
LIVERMORE,  CA  94550 

OICY  ATTN  DOC  CON  FOR  TECH  INFO  DEPT 
OICY  ATTN  DOC  CON  FOR  L-389  R.  OTT 
OICY  ATTN  DOC  CON  FOR  L-31  R.  HAGER 
OICY  ATTN  DOC  CON  FOR  L-46  F.  SEWARD 

LOS  ALAMOS  NATIONAL  LABORATORY 

P.O.  BOX  1663 

LOS  ALAMOS,  NM  87545 


OICY 

ATTN 

DOC 

CON 

FOR 

J. 

WOLCOTT 

OICY 

ATTN 

DOC 

CON 

FOR 

R.F.  TASCHEK 

OICY 

ATTN 

DOC 

CON 

FOR 

E. 

JONES 

OICY 

ATTN 

DOC 

CON 

FOR 

J. 

MALIK 

OICY 

ATTN 

DOC 

CON 

FOR 

R. 

JEFFRIES 

OICY 

ATTN 

DOC 

CON 

FOR 

J. 

ZINN 

OICY 

ATTN 

DOC 

CON 

FOR 

P. 

KEATON 

OICY 

ATTN 

DOC 

CON 

FOR 

D. 

WESTERVELT 

OICY 

ATTN 

D.  SAPPENFIELD 

SANDIA 

LABORATORIES 

P.O.  BOX  5800 

ALBUQUERQUE, 

NM 

87115 

OICY 

ATTN 

DOC 

CON 

FOR 

W. 

BROWN 

OICY 

ATTN 

DOC 

CON 

FOR 

A. 

THORNBROUGI 

OICY 

ATTN 

DOC 

CON 

FOR 

T. 

WRIGHT 

OICY 

ATTN 

DOC 

CON 

FOR 

D. 

DAHLGREN 

OICY  ATTN  DOC  CON  FOR  3141 

OICY  ATTN  DOC  CON  FOR  SPACE  PROJECT  DIV 

SANDIA  LABORATORIES 
LIVERMORE  LABORATORY 
P.O.  BOX  969 
LIVERMORE ,  CA  94550 

OICY  ATTN  DOC  CON  FOR  B.  MURPHEY 
OICY  ATTN  DOC  CON  FOR  T.  COOK 

OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
WASHINGTON,  D.C.  20545 

OICY  ATTN  DOC  CON  DR.  YO  SONG 


22 


OTHER  GOVERNMENT 

DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 
WASHINGTON,  D.C.  20234 
OICY  (ALL  CORRES: 

ATTN  SEC  OFFICER  FOR) 

INSTITUTE  FOR  TELECOM  SCIENCES 
NATIONAL  TELECOMMUNICATIONS  &  INFO  ADMIN 
BOULDER,  CO  80303 

OICY  ATTN  A.  JEAN  (UNCLASS  ONLY) 

OICY  ATTN  W.  UTLAUT 
OICY  ATTN  D.  CROMBIE 
OICY  ATTN  L.  BERRY 

NATIONAL  OCEANIC  &  ATMOSPHERIC  ADMIN 
ENVIRONMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COMMERCE 
BOULDER,  CO  80302 
OICY  ATTN  R.  GRUBB 
OICY  ATTN  AERONOMY  LAB  G.  REID 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 

AEROSPACE  CORPORATION 
P.O.  BOX  92957 
LOS  ANGELES,  CA  90009 
OICY  ATTN  I.  GARFUNKEL 
OICY  ATTN  T.  SALMI 
OICY  ATTN  V.  JOSEPHSON 
OICY  ATTN  S.  BOWER 
OICY  ATTN  D.  OLSEN 

ANALYTICAL  SYSTEMS  ENGINEERING  CORP 
5  OLD  CONCORD  ROAD 
BURLINGTON,  MA  01803 

OICY  ATTN  RADIO  SCIENCES 

AUSTIN  RESEARCH  ASSOC.,  INC. 

1901  RUTLAND  DRIVE 
AUSTIN,  TX  78758 
OICY  ATTN  L.  SLOAN 
OICY  ATTN  R.  THOMPSON 

BERKELEY  RESEARCH  ASSOCIATES,  INC. 
P.O.  BOX  983 
BERKELEY,  CA  94701 
OICY  ATTN  J.  WORKMAN 
OICY  ATTN  C.  PRETTIE 
OICY  ATTN  S.  BRECHT 


BOEING  COMPANY,  THE 
P.O.  BOX  3707 
SEATTLE,  WA  98124 
OICY  ATTN  G.  KEISTER 
OICY  ATTN  D.  MURRAY 
OICY  ATTN  G.  HALL 
OICY  ATTN  J.  KENNEY 

CHARLES  STARK  DRAPER  LABORATORY,  INC. 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02139 
OICY  ATTN  D.B.  COX 
OICY  ATTN  J.P.  GILMORE 

COMSAT  LABORATORIES 
LINTHICUM  ROAD 
CLARKSBURG,  MD  20734 
OICY  ATTN  G.  HYDE 

CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 
ITHACA,  NY  14850 

OICY  ATTN  D.T.  FARLEY,  JR. 

ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RICHARDSON,  TX  75080 
OICY  ATTN  H.  LOGSTON 
OICY  ATTN  SECURITY  (PAUL  PHILLIPS) 

EOS  TECHNOLOGIES,  INC. 

606  Wilshire  Blvd. 

Santa  Monica,  Calif  90401 
OICY  ATTN  C.B.  GABBARD 

ESL,  INC. 

495  JAVA  DRIVE 
SUNNYVALE,  CA  94086 
OICY  ATTN  J.  ROBERTS 
OICY  ATTN  JAMES  MARSHALL 

GENERAL  ELECTRIC  COMPANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GODDARD  BLVD  KING  OF  PRUSSIA 
P.O.  BOX  8555 
PHILADELPHIA,  PA  19101 

OICY  ATTN  M.H.  BORTNER  SPACE  SCI  LAB 

GENERAL  ELECTRIC  COMPANY 
P.O.  BOX  1122 
SYRACUSE,  NY  13201 
OICY  ATTN  F.  REIBERT 


GENERAL  ELECTRIC  TECH 
SERVICES  CO.,  INC. 

HMES 

COURT  STREET 
SYRACUSE,  NY  13201 
01CY  ATTN  G.  HILLMAN 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AK  99701 

(ALL  CLASS  ATTN:  SECURITY  OFFICER) 
01CY  ATTN  T.N.  DAVIS  (UNCLASS  ONLY) 
01CY  ATTN  TECHNICAL  LIBRARY 
01CY  ATTN  NEAL  BROWN  (UNCLASS  ONLY) 

GTE  SYLVANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A  STREET 
NEEDHAM,  MA  02194 

01CY  ATTN  DICK  STEINHOF 

HSS,  INC. 

2  ALFRED  CIRCLE 
BEDFORD,  MA  01730 

01CY  ATTN  DONALD  HANSEN 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
150  DAVENPORT  HOUSE 
CHAMPAIGN,  IL  61820 

(ALL  CORRES  ATTN  DAN  MCCLELLAND) 

01CY  ATTN  K.  YEH 

INSTITUTE  FOR  DEFENSE  ANALYSES 
1801  NO.  BEAUREGARD  STREET 
ALEXANDRIA,  VA  22311 
01CY  ATTN  J.M.  AEIN 
01CY  ATTN  ERNEST  BAUER 
01CY  ATTN  HANS  WOLFARD 
OlCY  ATTN  JOEL  BENGSTON 

INTL  TEL  &  TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

OlCY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

11011  TORREYANA  ROAD 
P.O.  BOX  85154 
SAN  DIEGO,  CA  92138 

OlCY  ATTN  J.L.  SPERLING 


JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 

OlCY  ATTN  DOCUMENT  LIBRARIAN 
OlCY  ATTN  THOMAS  POTEMRA 
OlCY  ATTN  JOHN  DASSOULAS 

KAMAN  SCIENCES  CORP 
P.O.  BOX  7463 

COLORADO  SPRINGS,  CO  80933 
OlCY  ATTN  T.  MEAGHER 

KAMAN  TEMPO-CENTER  FOR  ADVANCED  STUDIES 
816  STATE  STREET  (P.O  DRAWER  QQ) 

SANTA  BARBARA,  CA  93102 
OlCY  ATTN  DASIAC 
OlCY  ATTN  WARREN  S.  KNAPP 
OlCY  ATTN  WILLIAM  MCNAMARA 
OlCY  ATTN  B.  GAMBILL 

LINKABIT  CORP 
10453  ROSELLE 
SAN  DIEGO,  CA  92121 
OlCY  ATTN  IRWIN  JACOBS 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC 
P.O.  BOX  504 
SUNNYVALE,  CA  94088 
OlCY  ATTN  DEPT  60-12 
OlCY  ATTN  D.R.  CHURCHILL 

LOCKHEED  MISSILES  &  SPACE  CO.,  INC. 

3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

OlCY  ATTN  MARTIN  WALT  DEPT  52-12 
OlCY  ATTN  W.L.  IMHOF  DEPT  52-12 
OlCY  ATTN  RICHARD  G.  JOHNSON 
DEPT  52-12 

OlCY  ATTN  J.B.  CLADIS  DEPT  52-12 

MARTIN  MARIETTA  CORP 
ORLANDO  DIVISION 
P.O.  BOX  5837 
ORLANDO,  FL  32805 
OlCY  ATTN  R.  HEFFNER 

M.I.T.  LINCOLN  LABORATORY 
P.O.  BOX  73 
LEXINGTON,  MA  02173 

OlCY  ATTN  DAVID  M.  TOWLE 
OlCY  ATTN  L.  LOUGHLIN 
OlCY  ATTN  D.  CLARK 
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MCDONNEL  DOUGLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
01CY  ATTN  N.  HARRIS 
01CY  ATTN  J.  MOULE 
01CY  ATTN  GEORGE  MROZ 
01CY  ATTN  W.  OLSON 
01CY  ATTN  R.W.  HALPRIN 
01CY  ATTN  TECHNICAL  LIBRARY  SERVICES 

MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
01CY  ATTN  P.  FISCHER 
01CY  ATTN  W.F.  CREVIER 
01CY  ATTN  STEVEN  L.  GUTSCHE 
01CY  ATTN  R.  BOGUSCH 
01CY  ATTN  R.  HENDRICK 
01CY  ATTN  RALPH  KILB 
01CY  ATTN  DAVE  SOWLE 
01CY  ATTN  F.  FAJEN 
01CY  ATTN  M.  SCHEIBE 
01CY  ATTN  CONRAD  L.  LONGMIRE 
01CY  ATTN  B.  WHITE 

MISSION  RESEARCH  CORP. 

1400  SAN  MATEO  BLVD.  SE 
SUITE  A 

ALBUQUERQUE,  NEW  MEXICO  87108 
01CY  R.  STELLINGWERF 
01CY  M.  ALME 
01CY  L.  WRIGHT 

MITRE  CORPORATION,  THE 
P.O.  BOX  208 
BEDFORD,  MA  01730 

01CY  ATTN  JOHN  MORGANSTERN 
01CY  ATTN  G.  HARDING 
01CY  ATTN  C.E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 
01CY  ATTN  W.  HALL 
01CY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
12340  SANTA  MONICA  BLVD. 

LOS  ANGELES,  CA  90025 

01CY  ATTN  E.C.  FIELD,  JR. 


PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(NO  CLASS  TO  THIS  ADDRESS) 

01CY  ATTN  IONOSPHERIC  RESEARCH  LAB 


PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  3027 
BELLEVUE,  WA  98009 

01CY  ATTN  E.J.  FREMOUW 

PHYSICAL  DYNAMICS,  INC. 

P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 

R  &  D  ASSOCIATES 
P.O.  BOX  9695 
MARINA  DEL  REY,  CA  90291 
01CY  ATTN  FORREST  GILMORE 
01CY  ATTN  WILLIAM  B.  WRIGHT,  JR 
01CY  ATTN  ROBERT  F.  LELEVIER 
01CY  ATTN  WILLIAM  J.  KARZAS 
01CY  ATTN  H.  ORY 
01CY  ATTN  C.  MACDONALD 
01CY  ATTN  R.  TURCO 
01CY  ATTN  L.  DeRADD 


RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
01CY  ATTN  CULLEN  CRAIN 
01CY  ATTN  ED  BEDROZIAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 
01CY  ATTN  BARBARA  ADAMS 

RIVERSIDE  RESEARCH  INSTITUTE 
80  WEST  END  AVENUE 
NEW  YORK,  NY  10023 

01CY  ATTN  VINCE  TRAPANI 


PHOTOMETRICS,  INC. 

4  ARROW  DRIVE 
WOBURN,  MA  01801 

01CY  ATTN  IRVING  L.  KOFSKY 
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SCIENCE  APPLICATIONS,  INC. 

P.O.  BOX  2351 
LA  JOLLA,  CA  92038 

01CY  ATTN  LEWIS  M.  LINSON 
01CY  ATTN  DANIEL  A.  HAMLIN 
01CY  ATTN  E.  FRIEMAN 
01CY  ATTN  E.A.  STRAKER 
01CY  ATTN  CURTIS  A.  SMITH 
01CY  ATTN  JACK  MCDOUGALL 

SCIENCE  APPLICATIONS,  INC 
1710  GOODRIDGE  DR. 

MCLEAN,  VA  22102 
ATTN:  J.  COCKAYNE 

SRI  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 

01CY  ATTN  DONALD  NEILSON 
01CY  ATTN  ALAN  BURNS 
01CY  ATTN  G.  SMITH 
01CY  ATTN  R.  TSUNODA 
01CY  ATTN  DAVID  A.  JOHNSON 
01CY  ATTN  WALTER  G.  CHESNUT 
01CY  ATTN  CHARLES  L.  RINO 
OiCY  ATTN  WALTER  JAYE 
01CY  ATTN  J.  VICKREY 
OICY  ATTN  RAY  L.  LEADABRAND 
OICY  ATTN  G.  CARPENTER 
OICY  ATTN  G.  PRICE 
OICY  ATTN  R.  LIVINGSTON 
OICY  ATTN  V.  GONZALES 
OICY  ATTN  D.  MCDANIEL 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  MA  01730 

OICY  ATTN  W.P.  BOQUIST 

TOYON  RESEARCH  CO. 

P.O.  Box  6890 
SANTA  BARBARA,  CA  93111 
OICY  ATTN  JOHN  ISE,  JR. 

OICY  ATTN  JOEL  GARBARINO 


TRW  DEFENSE  &  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90278 
OICY  ATTN  R.  K.  PLEBUCH 
OICY  ATTN  S.  ALTSCHULER 
OICY  ATTN  D.  DEE 
OICY  ATTN  D/  STOCKWELL 
SNTF/1575 

VISIDYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MASS  01803 
OICY  ATTN  W.  REIDY 
OICY  ATTN  J.  CARPENTER 
OICY  ATTN  C.  HUMPHREY 


DOE  DISTRIBUTION 


University  of  California 
Lawrence  Livermore  National  Lab 
Post  Office  Box  808 
Livermore ,  CA  945S0 
H.G.  Ahlstrom,  L-481 
J .L  Emmett ,  L-448 
J.F.  Holzrlchter,  L-481 
M.J.  Honsler,  L-479 
J.H.  Nuckolls,  L-477 
L.W.  Coleman,  L-473 
J.T.  Hunt,  L-481 
A.B.  Langdon,  L-477 

U.S.  Department  of  Energy 
Office  of  Inertial  Fusion 
Washington,  DC  20545 
L.E.  Killian 
G.  Gibbs 
T.F.  Godlove 

S. L.  Kahalas 
J.E.  Lewis 

R. L.  Schriever 

T. H.  Walsh 

S. J.  Barish 

U.S.  Dept,  of  Energy  (194  cys) 
Technical  Information  Center 
P.O.  Box  62 
Oak  Ridge,  TN  37830 

Defense  Tech.  Information  Ctr.  (2  cys) 
Cameron  Station 
5010  Duke  Street 
Alexandria,  VA  22314 

R.  McCrory 

University  of  Rochester 
250  East  River  Road 
Rochester,  NY  14623 

Robert  T.  Duff 

U.S.  Department  of  Energy 
Office  of  Classification 
Washington,  DC  20545 


Rex  B.  Purcell  (2  cys) 
U.8.  Department  of  Energy 
Nevada  Operations  Office 
Post  Office  Box  14100 
Us  Vegas,  NV  89114 


LIST  FOR  REPORTS 


Z.N.  Zafiris/R.  Bredderman 
U.S.  Department  of  Energy 
San  Francsico  Operations  Office 
1333  Broadway 
Oakland,  CA  94512 

Los  Alamos  National  Uboratory 
Post  Office  Box  1663 
Los  Alamos,  NM  87545 

S.D.  Rockwood,  ICF  Prog.  Mgr 
DAD/IF,  M/S527  (6  cys) 

G.  Yonas  (4  cys) 

Sandia  National  Uboratorles 
Post  Office  Box  5880 
Albuquerque,  NM  87185 

S .  Bodner 

Naval  Research  Laboratory 
Code  4730 

Washington,  DC  20375 

T.  Coffey 

Naval  Research  Laboratory 
Code  1001 

Washington,  DC  20375 

Alexander  Glass 
XMS  Fusion,  Inc. 

3941  Research  Park  Drive 

P.O.  Box  1567 

Ann  Arbor,  MI  48106 

NRL  Code  4700  (26  cys) 

NRL  Code  4730  (200  cys) 

NRL  Code  2628  (22  cys) 


27 


